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1. 


Introduction 


Suppose  a  random  sample  of  points  P^,  i=l,  2,  ,  N  ,  is  given 

on  the  surface  of  a  sphere  with  center  0  and  radius  1  .  The  vectors 
OP.  may  denote  directions.  In. some  geophysical  problems  vectorial  data 
of  this  type  occur,  in  .which.. the.  vectors ..OP  are  more  or  less  coplanar, 
that  is,  clustered  close  to  a  plane  through  the  center  0  .  Such  a  plane 
intersects  the  sphere  in  . a  great.circle,  which  may  be  regarded  as  the 
equator  of  the  sphere;  the  probability: distribution  over  the  Spherical 
surface  which  is  then  appropriate  .  to  describe  the  distribution  of  points 
P  (or  equivalently  of  unit. vectors  OP)  must  have  high  density  near  the 
equator  and  low  density  at  the  North  and  South  Poles,  Such  a  distribution 
is  called  an  equatorial,  distribution . - ~  Suppose  spherical  polar  coordinates 
(0,  <£)  are  used  for  P  .,  and  let  0  =  0  be  the  North  Pole;  here 
0  <■  4>  _<  2?r  and  0  <_  0  IT  .... The.. density 

2 

(1)  f(6,  <p).=  c3  (k)  sin  0  exp  (-Kf  cos  0)  , 

with  K  >  0  and  c^(k)  a  normalizing, constant,  was  introduced  as  a 
suitable  equatorial  distribution::by  Watson  (1965),  If  K  <  0  ,  the  den¬ 
sity  may  be  used  to  describe, a: bimodal  population,  with  modes  in  the 
density  at  the  poles  and  low  density, on  the  equator  0  =  ir/2  .  Clearly, 
suitable  orthogonal  transformations  will  place  the  polar  axis  along  any 
vector  A . «.  When  K  =  0  ,  .the  .density  reduces  to  the  uniform  density 
over  the  sphere, 

A  common  problem,  given, sample, data,  is  first  to  test  for  uniformity, 
that  is,  to  test  the  null, hypothesis  H^:  K  =  0  against  the  alterna¬ 
tive  H^:  K  >  0  or  H^i.-K  <0,,  For  the  equatorial  alternative, 

,  Watson  (1965)  showed  that  .the  likelihood  ratio  test  of  Hq  , 
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when  the  antimodal  axis  A  is  known,  is  based  on  S  =  N  ^  5  cos^  a  , 

~  * -  i=l  i 

where  a_^  is  the  angle  between  •  ~A  and  0P_^,  i=l,  2,  ...  ,  N  .  The 
vector  A  may  point  to  either,  antimode .  When  S  is  significant  in  the 
lower  tail,  Hq  is  rejected.. in.. favor  of  ...  S  may  also  be  used  to 

test  for  uniformity  against.. the  bimodal.  density  given  by  (1)  with  K  <  0  ; 
then  Hq  is  rejected  in  favor  of  if  S  is  significant  in  the  upper 

tail.  Percentage  points  of  -S. -(there  called  T)  are  given  in  Stephens 
(1966). 

In  this  paper  we  give  significance  points  for  the  test  statistic  to 
be  used  for  testing  Hq  against  H^  and  for  testing  Hq  against  Hg 
in  the  more  realistic  case  when  the  antimodal  and  modal  vectors  A  , 
respectively,  are. not  known.  We  also  give  some  indications  of  the  power 
of  the  test.  The  problems  posed  in  finding  the  distributions  of  the  test 
statistics  have  interesting  connections  with  principal  component  analysis 
and  with  an  engineering  application;  they,  have  not  been  solved  exactly 
except  in  one  special  case, -but-signif icance  points  are  given,  based  on  a 
combination  of  asymptotic  results  and  Monte  Carlo  studies.  We  conclude 
the  paper  with  an  explanation  of .  the-corresponding  problem  in  two  dimen¬ 
sions,  because  there  the  distribution  theory  may  be  solved  exactly. 

2.  Tests  for  Randomness  Against - Equatorial  or  Bimodal  Alternatives 

2.1.  The  procedures.  We  now  discuss  the  test  of  Hq  when  the- vector 
A  ,  pointing  to  an  antimode  for.  the  equatorial  alternative  and  to  a  mode  for 
the  bimodal  alternative,  is  not  known.  The  likelihood  ratio  test  statistic 
is  found  as  follows  (Watson,  1965).  Suppose.for  a  given  vector  A'  , 

S'  =  N  )  .  .  cos  a.  ,  where  a.  is  the  angle  between  OP.  and  A'  . 

-  1=1  1  1  ^ 
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As  A'  is  allowed  to  vary. over  the  surface  of  the  sphere,  S'  will  take 

values  between  a.minimum,  S  „  ,  and  a  maximum,  S  .  S  ,  will  be 

min  max  mm 

the  likelihood  ratio  test  statistic  for  the  test  of  H„  against  H 

U  ii 

(equatorial  alternative)  and  S  will  be  used  against  H,,  (bimodal 

max  & 

alternative).  S  .  and  S  may  easily  be  calculated  from  the  latent 
min  max 

roots  (eigenvalues,  characteristic  roots)  of  a  matrix  Q  to  be  given 
below.  The  tests  may  therefore  be-set. out  as  follows: 

(a)  Suppose  £_^,  nu,  n_^  are  the  direction  cosines  of  OP^,  referred 

to  a  fixed  system  of  coordinates;  that  is,  let  =  sin  0_^  cos  (|)^  , 

m.  =  sin  0.  sin  <t> .  ,  n.  =  cos  0.  ,  i=l,...„  ,  N  . 

1  l  i  i  i 

(b)  Let  Q  be  the  following  matrix  divided  by  N  : 


N 


) _  a: 

i=i  1 

N 

5  A.  m. 

(-  •  ■  -i  -i 

i=l 


N 


i=l 


h  n± 


N 


5  A,,  m. 

. - T-  1  1 

1-1 


N 


i=l 

N 


i=l 


m , 


mi  ni 


N 


YZ  K  n- 

.~r  i  i 
1=1 


N 


i=l 


N 


n. 


i=l 


(c)  Let  the  latent  roots  of  Q  be  .  _>  >_  ,  and  let 


the 


corresponding  latent  vectors  be  u, ,  u„,  and  u„  Then  S  ,  =  A„  and 

~1  ~2  ~3  mm  3 

S  =  A,  . 
max  1 

(d)  If  is  less  than. the  appropriate  entry  for  values  of  N 

and  a  in  Table  1,  reject  Hq  at  significance  level  a  in  favor  of  the 
equatorial  alternative;  in  this  event  u^  is  the  maximum  likelihood  esti¬ 


mate  of  the  antimodal  vector  of  the  distribution.  Alternatively,  if  S 


max 


is  greater  than  its  appropriate. entry,  reject  Hq  in  favor  of  the  bimodal 


alternative,  with  u^  as  the  maximum  likelihood  estimate  of  the  modal  vector. 
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2.2  Related  problems. 

(a)  Principal  component. analysis.  Let  the  vectors  OP^  be  denoted 


=  2  .  = 

2,  1 

= 

2i 

hi  i 

/  ^ 


l  > 


n. 


\  X 


i=l,  ...  ,  N  . 


A  vector  3  =  (3^>  &2’  £3)  *  normalized  s0 


(2) 


3’3  = 


k=l 


gk=  1 


has  direction  cosines  as  components.  Then 


r  =  b'  z  =  3 

1  ~  ~1  , — r- 


k=l 


k  zki 


is  the  signed  length  of  the  projection  of  OP  on  3  ;  in  other  words  it 
is  the  length  of  OP^  in  the  direction  of  3  .  The  sum  of  squares  of 


these  lengths  is 
N 


N 


(3)  L_  vj  - 

i-1  i=l 


O  N 

(3'  z.r  =  y~  3'  Z.  z!  3  =  N3VQ3 

~  ~  i  ,  i  ~  ~  i  ~  l  ~  ~  ~  ~ 

l-l 


The  maximum  of  (3)  with  respect  to  3  ,  subject  to  (2),  is  the  largest 

Thus  N  S 


characteristic  root  of  NO  ,  which  is  NX.,  =  N  S 

~  1  max 


max 


is 


the  sum  of  squared  distances  of  the  points  P^,  ...  ,  P^  maximized  with 
respect  to  direction.  The  maximizing  -  direction  3^^  defines  the  first 
(sample)  principal  component  (from  the  origin)  as  v^^  =  3^^  z  ,  and 

s  ,-  i  e:  v«2 

max  N  .  v"  i 
i=l 

is  the  (sample)  variance  (about  0)  of  the  first  principal  component. 
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The  projection  of  OP.  -  on.  the  plane  . or thogonal  to  3^^  is 

~~  1  ~ 

z  -  8^"^  •  The.  signed  length -of  the  projection  of  this  vector  on 

(1) 

a  normalized  vector  3  in  this -plane  is  3’  z.  because  3'3V  =  0  . 

rsj  A/  1  /v 

The  sum  of  squares  is  (3)  and . the  maximum  is  the  second  characteristic 

(2) 

root  of  NQ  ,  which  is  NX„  ...  ..Let  the  maximizing  direction  define  3 

■V  4  ~ 

and  let  8^^  be  orthogonal  to  8^^  and  8^^  .  Then  Q3^^ 

(3) 

is  the  sum  of  squares  of  lengths  in  the  direction  of  3  ,  and  this  is 

the  minimum,  namely  NA^  =  N  . 

The  discussion  here  parallels  the  treatment  of  principal  components 
in  Section  10.2  of  Anderson  (1958),  with  expectation  replaced  by  sample 
average  (or  sum) ,  coordinates  from  the  mean  replaced  by  coordinates  from 
the  origin,  and  p  set  at  3  . 

The  density  (1)  can  be .. obtained  from  a  three-dimensional  normal 
distribution  with  means  0  ,  covariances  0  ,  and  two  variances  equal. 

Then  (1)  is  the  conditional  density  of  the  two  angles  in  polar  coordinates 
for  the  radius  fixed. 

(b)  Mechanical  Engineering. --If -unit  masses  are  located  at  each 

point  P^  ,  the  quantity  N(l-S)  gives -the  moment  of  inertia  about  A'  , 

and  N(l-S  .  )  and  N(l-S  )  are  -respectively  the  maximum  and  minimum 
mm  max  r 

moments  of  inertia  obtainable  by  varying  Af  .  Thus  the  distributions 

of  S  .  and  S  can  be  used  to  give  the  distributions  of  maximum 
mm  max 

and  minimum  moment  of. inertia  obtainable  from  a  random  distribution  of 
N  unit  masses  on  the. sphere. 
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3. 


Distributions  of  ST  . -and ..of  S 

_  - . . mm  — : .  -  -  max 

3.1.  A  large-sample - result ...  The  exact  distributions  of  S  .  and 

mm 

of  are  very  difficult  to. find. :  In  this  section  we  give  large- 

sample  results,  leading  to  very. good -approximations  to  significance  points, 

-h 

to  the  order,  of  N  .  As  N..-K®0  >  ..the  entries  in  Q  have  an  asymptotic 

joint  normal  distribution.  From  it. we  can  obtain  the  asymptotic  joint 

distribution  of  the  latent,  roots;  ..the  marginal  distributions  of  the  largest 

and  smallest  roots  are  the  asymptotic  distributions  of  S  and  S  . 

max  mm 

Since  we  want  the  asymptotic  distribution  for  the  two-dimensional  case  as 
well  as  the  three-dimensional  case  and  it  may  be  useful  for  higher  dimen¬ 
sional  cases,  we  prove  a  more  general -result. 

Suppose  y  =  (y^  ...  ,  y^)  '  has  a  p-dimensional  normal  distribution 
with  mean  0  and  covariance  matrix  I  .  Let  z  =  y//y’y  .  Then  z  has 

the  uniform  spherical  distribution. in .  p  dimensions,  and  z  and  w  =  /y 'y 

2  2 

are  independent.  Since  w  has  a  X  -distribution  with  p  degrees  of 
freedom, 


*>  — m  *  ,.2slga  „%m  r[Jg(m+p)] 

-  £(Xp)  -2  -fgf- 


Then  0  =  |>y  =  ^Lwz  =  ^w  implies  z  =  0  . 
independent,  a  moment  of  the  components  of  z  is 
of  the  components  of  y  divided  by  the  moment  of 
Let  A  =  (a^)  =  zz'  ;  then 


Because  z  and  w  are 
the  corresponding  moment 
w  of  the  same  degree. 


%A  =  £yy’  -  ~  I  , 

~  %Mt.  P  ~ 

2  _  ^  4  =  1  ^  4  .  .  1 

akk  4~k  ,.2  ^yk  p(p+2)  ’ 

$w  ~  r  r 


6 


1 


k  £  ft  , 


~e  2  _  ^  ?  2  2 

^  ak£  ~  ^aivaoV  ~  ^zvzc 


lkk  ftft 


Jk  ft 


^  akkak£  5  zkz£ 


U  w 
1 

-e  ^ 


2.2  1 

4  ?ykyft  =  p  (p+2)  * 


£ykyft  =  0  ’ 


The  ..variances -and  ~covariances» -of -the  .elements  of  '  A  are 


Var  akk  " 


to:Vu> 


P (p+2)  ^2 


•=  2  - 


-P--1 


P (p+2)  p2 


P  (p+2) 
2 

p2(p+2) 


Var  <akJi) 


P (p+2)  5 


k  ^  A  . 


k  t5  A  , 

k  ^  &  , 


and  0  for  other  covariances. 

Let  z.,  i“l,  2,  . ...  ,  N,  be  the  sample  values  of  z  ,  and  define 

~  1  ~ 

N 

Q  =  (1/N)  YZ  z.zi  > 

i«l  ~1~1 

X  =  /N  (Q  -  -  I)  . 

~  ~  P  ~ 

T -  N 

The  trace  of  Q  is  (1/N)  /  .  ..  z!z.  =  1  ,  so  the  trace  of  X  is  zero. 

~  * 1=1  ~1~1  ~ 

Then  X  has  a  limiting  normal  distribution  with  mean  0  and  the 
same  covariances  as  A  .  In  *  the  limiting  normal  distribution  of  X  , 
the  (functionally '.independent)— off-diagonal  terms  are  mutually  (statis¬ 
tically)  independent -and -independent  of  the  diagonal  terms.  Let  v’ 

Jc  -  (x. ,  x00,  ...  ,  x  .  )  1  be  the  vector  of  diagonal  terms;  its  covari- 
~  -  11  zz  pp 

ance  matrix  is 
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2 

p2(p+2) 


r: 


-i 

p-i 


-i 


-i 

-i 

« 

p-i 


2 

p2(p+2) 


[pl-ee*]  , 


where  £  =  (1,  1,  ...  ,  1)'.  This  is  a  singular  matrix  and  the  limiting 

distribution  of  X  is  singular..  In  order  to  deal  with  densities,  we 

study  first  the  latent  roots  of  the  matrix  II  =  X  +  vl  ,  where  v  is 

2 

normally  distributed  with  mean  0  and  variance  2/ [p  (p+2)]  .  From  the 
density  of  the  latent. roots  of  U  ,  we  can  obtain  the  density  of  any 
p-1  roots  of  X  . 

Let  u  =  (u. .  ,  u,.,. ,  .  .  .  ,  u  . ) *  =  x  +  v£  .:  With,,  v  chbsett  With  the 

~  JL1  Z.L  PP  ^  ~ 

above  variance,  the  covariance  matrix  of  U  is  diagonal: 


%  uu*  =  - -  [pl-es’]  +  -r  1 


ee 


T  - 


P  (p+2) 


p  (p+2) 


P  (p+2  )  ~ 


I  . 


The  variance  of  u^  is  2/ [p (p+2) ]  ,  which  is  twice  the  variance  of 
u^,  k  ^  Z  ;  the  latter  is  the  variance  of  a^  ,  given  above.  The 
density  of  this  limiting  distribution-of  U  is  then  proportional  to 


'■ 


t 

(  rh 

fl  -4  n  ^ 

k=l  .  k<Z 

\ 

/  1 

i.  P (p+2)  V2-  2 

2  +  1 
p(p+2)  p (p+2 ) 

*  2  + —  uU 

Z  k,A-l 

J 

k  i 

[p(p+2)]p<p+1)/4  -Js[p(p+2)/2]tr  D2  . 

2lip  (2p)p(p+1)M  "j  ~ 

If  the  roots  of  U  with  the  above  density  are  s^  >  s£  >  ...  >  s^  , 

their  density  is  [by .Theorem  .13.3.1  of  Anderson  (1958)] 
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[ P (p+2 ) ] 


p (p+1/4 


f 


2P(p+3)/4  jEj-  r[Mp_1+1)]  [ 

i=l 


exp  <J-%[p(p+2)/2]  sH  TJ  (s  -s  )  , 

1=1  1  i,j=l  J 


i<j 


s  ^  •  •  •  <  S-  • 
p  -  -  1 


Now  the  roots  r.  >  r„  >  ...  >  r  of.  X  are  related'!t6  those  of  U 
1  -  2  '  P  ~ 

(S1  >  s2  >  >  Sp)  by  the-  relation 


x.  =  r.  +  v  . 

l  l 


Since  tr  X  =  0  ,  we  have  \  p  r.  =0  ;  then  pv  =  3  p  s.  . 

~  £ -  i=l  i  * -  i=l  i 


So  v=  (s^  +  . ..  +s  )/p  and 


P-1  1  7“ 

ri  =  si-v  =  Vsi-F  W  v 


1—1 f  • • •  y  p— 1  • 


j=l 


Then 


± 

•  _  1  x  •  _  1 


2  ,  2 
ri  +  pv  , 


i=l  *  i=l 

where  r  =  -(r^  +  ...  +  r  _^)  .  The  joint  density  of  v  and 


r1  ,  ...  ,  r  is 
1  p-i 


jP  (p+1 )  /4  (p+2  )  p  (p+1)  lk-h  f 

,p(p+3)/4- 


(4) 


1  i-1  71  J 


X  e-Jsl»2<P+2)/21v2-  ,  r  .  -r.  -  -  r  ,  .<  r  .<...<  r. 

2  7T5  "  ’pl  -  P  -1  ^  P  -1  1 


Thus  v  is  independent  of  ...  ,  r  ,  and  the  density  of  the  latter 

set  is  the  first  factor.  The  density  of  the  limiting  distribution  of  the 

p-1  largest  characteristic  roots  of  v'NCQ  -  —  I)  is  the  first  factor 

P  ~ 

and  the  density  of  the  limiting  marginal  distribution  of  any  subset  of 
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these  characteristic  roots- is -obtained .as -a  marginal  density  from  the 
factor .  [ Justification. follows  from. the. Theorem  on  Limiting  Distribution, 

p.  140,  Anderson  (1963).] 

The  three-dimensional- case.  When  p=3  ,  the  above  gives  the  density 
of.  the  limiting  distribution  of-.r^.  and  r^  »  Our  interest  is  in  r^ 
and  r^  =  -r^  -x^  ■>  and  the  density  of  the  limiting  distribution  of  these 
variables  is 


152  /45 
4  /2t r 


exp 


15 


2  2 

r  4*  r  4*  r  r 
13  1  3 


(2rL  +  r3)(r1  -  r3) (-r1  -  2r3) 


-2r3  £  r3  <_  ~hx-^  £  0 


Integration  with  respect  to  r3-  gives  the  density  of  r^  as 

2  2 
/45 


(5)  f (r  )  = 

/27T 


135  2  -arl  -arl  ^  -4arl  1 

r,  e  -e  +  e 


0  <  r  < 


8  1 

where  a  =  45/8  .  The  density  of  r3  is  the  same  but  for  r3  _<  0  .  The 
limiting  distribution  of  x^  ,  found  from  above,  is  normal  with  mean  0 
and  variance  1/45  .  ' 

For  calculations,  it  is  convenient  to  write  the  distribution  of  r^ 
in  terms  of  u  =  (r^  /A5)/2  ;  then  integration  of  (5)  gives 


Pr 


(u  >  z)  =  Pr(r1  >  2z//45)  =  — ze  z  +  2  -  $(z)  -  $(2z)  , 
1  /2tt 


2  2 

where  $(z)  =  -  /  e  r  dt  ,  and  z  >  0  .  Similarly, 


2 

Pr(r„  <  2z/v/45")  =  (-ze  Z  /2 

3  /2i 


+  $(z)  +  $(2z)  , 


z  <  0 


The  critical  values*  of . .r^  .,  .at .  significance  level  a  ,  are  (r^)  : 


a: 

.10  .  .08 

_  .05  . 

.025- 

.01 

.005  .002 

<rl> a 1 

0.788  0.816 

0.873  - . . 

0.948 

1.038 

1.100  1.177 

Critical 

values  of  r^ 

are  (r_)  = 

.  -  -  3  a 

.  Then 

critical  values 

(S  .  ) 
mm  a 

and  (S  )  are,  to  order 

max  a 

N~^  : 

(6) 

(  Wc.  =  I  + 

i 

-H|  CO 

II 

s ; 

>  H: 

S-/ 

<rl>a  . 

(S  ) 
max  a 

-  1  . 

3 

:3«2.  :  A  sitall-sample  .result.  .  Let^jfZ  =  '(z,  .  >  z  ) ,  where  'z:  iis 

defined  above.  Then  .  NQ  .=  ZZ'.'.:and  the1  latent’ roots  of  ZZ'  are  the  latent 
roots  of  Z'Z"  except:for  ,0.-r66ts*  equal  ..in  number  to  the  difference ’iri 
.dimensionalities ;.of  ZZ'  and.  Z'Z;  .:  Then:- the' (nonzero)  latent  roots  of 


NQ 


are  the  (nonzero)  latent  toots  of 


Ztli 


say.  R  is  an  N  x  N  correlation  matrix  based  on  p  observations  from 
N(0,  I)  ,  using  deviations  from  -0-,  instead  of  the  mean. 

For  the  particular  case  p=3  ,  if  N=2  or  3  ,  R  has  a  density;  if 
N  >  3  ,  R  does  not  have  a  .density  because  its  rank  is  3  which  is  less 
than  its  dimensionality  - N. . . . If  N=2 


*  We  are  indebted  to  R.L.  Anderson  for  carrying  out  these. computations. 
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and  under  the  null  hypothesis  r  has  the  uniform  density  over  (-1,  1)  . 

If  FN3  , .  the  density  of  ..r^*  rl3*  r23  ^S'  ^  for  R  positive 

definite  and  0  otherwise* . See . Anderson  (1958),  pp.  64  and  175.  The 

number  of  degrees  of  -  freedom. .  n  is  3  -  because  here  deviations  from  0 

are  used,  instead  of  from . the.  mean ; . the  dimensionality  p  is  N  .] 

When  N=2  ,  the . characteristic  roots  of  R  are  1+r  and  1-r  . 

Thus  S  =  %[1  +  max(r,  -r]  and  S  ,  =0  because  Q.  is  singular, 

max  .  mm  ~ 

The  cdf  of  S  is  (for  %  <  x  <  1) 
max  —  — 

Pris  <  x}  .=  Pr{max(r ,  ~r)  <  2x  -  1}  =  2x  -  1 
max  —  — 

and  the  significance  level  based. on. a  test  with  significance  point  x  is 
2(1  -  x)  .  Some  values  are  given -below. 


Table  3.1 


Probability 

.90  . 

.95 

.975 

.99 

Significance  level 

.1 

.05 

.025 

.01 

Significance  point 

.95 

.975 

.9875 

.995 

3.3.  Significance  Points.  For  smaller  values  of  N  ,  exact  critical 

values  of  S  .  and  S  will  increasingly  depart  from  the  values  in 

(6).  These  exact  values  . have  been  estimated  by  Monte  Carlo  studies,  using 

mostly  10,000  samples ,. each  of  size  .N  ,  . for  N  =  4(1)12,  14,  15,  16,  20, 

30,  50,  60,  80  for  10,000.  samples,  for  N  =  5,  10,  20,  50  for 

S  and  the  middle  root  A0_ ...  The  Monte  Carlo  points,  for  given  a  were 
iricix 
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— *< 

plotted  against  N  ;  both  the. lower  tail  points  of  ,  and  the 

upper  tail  points  of  S  approach- smoothly,  from  above,  the  asymptotic 

mcix 

(6) *  Smoothed  curves  through . these,  points . give  the  points  in  Table  1. 

Analysis  of  these. curves  suggests  that  these  percentage  points  are  accurate 

to  within  one  or  two  units -in  the -third  decimal  place,,  For  N  beyond 

the  table  the  expressions  in- (6)  give  a  very  good  approximation.  The 

critical  values  of  S-.  and  S.  ,  for  small  N  ,  are  not  symmetric 

with  respect  to  1/3  as  are  the .asymptotic  values;  this  is  to  be  expected 

as  S  .  has  range  0  to  1/3- ,  while  S  has  range  1/3  to  1  .  For 
mm  max  ° 

the  middle  latent  root  of  Q  ,  Monte  Carlo  values  agree  asymptotically 

with  0.3333  -  (r^)  /v^N-  ,  where,  r^  is  distributed  N(0,  1/45)  ;  for  smaller 
N  ,  the  distribution  becomes .gradually  unsymmetric.  Some  critical  values 
are  given  in  Table  2. 

3.4.  Power  Studies.  Table  3  gives  an  indication  of  the  power  of  the 
test  for -randomness,  when  applied -at  level  a  ,  and  when  the  sample  is  of 
size  N  and  in  fact  drawn . from.an  equatorial  distribution  with  a  given 
K  .  Thus,  for  example,  when. the  test -is  applied  at  the  5%  level,  with 
N  =  20  and  K  =2  ,  the  power  - (i.e„ , .probability  of  declaring  a  sample 
significant,  i.e„,  not- random) -is  0.400  ;  if  the  sample  is  of  size  10 
and  K  =  3  ,  the  power  is  0.300  . 

Monte  Carlo  results  were  used  to  find  the  proportion  of  1,000  samples 
(for  N  <  20)  or  400  samples  (for  N  =?  20,  30,  60)  which  were  significant 
when  the  test  was  applied. as -described.-  Curves  of  proportion  significant 
were  plotted  against  l/N.,_for  given.  K  and  a  ,  and  used  to  construct 
the  table. 


13 


4. 


The  Two-Dimens ional - Prob lem 


It  is  interesting -to  observe -that -the  .distributions  of  S  .  and  of 

mm 

S  can  be  solved_exactly-when-the.points  P.  are  on  the  circumference 
max  - -  1 

of  a  circle.- -Suppose. .  Q^-.-dLs -the.- polar  coordinate  of  P/  referred  to  a 
suitable  origin,  and. let  the radius  be  1  ;  in  its  most  general  form,  the 
distribution-corresponding  to . (1)  is 

(7)  f  (6 )  =  c2(<)  exp  [-K.  cos2  (0--  0Q)J  ,  0  £  0  _<  2tt  . 


This  distribution .has  two  modes  and  two  antimodes  symmetrically 
placed  round- the  circle;  when  K  <_  0  ,  0  -  0^  and  0  =  0^  +  TT  gives  the 
antimodes,  and  when  K  >  0 ... ,  they .  give  the  modes .  A  change  in  sign  of 
K  is  exactly  equivalent -to-a  rotation  of  coordinates  through  90°  ,  so 
that  no  distinction -can -be  made  between. the  equatorial  and  bimodal  cases, 
and  a  test  of  uniformity  against  -  the  alternative  (7)  should  give  the  same 
result  whether,  for  the  the-alternative,  K  is  positive  or  negative.  We 
can  see  that  this  is  so -as  follows. ...  -  The  test  statistic,  when  A  is  not 
known ,  is 


mm 


1  . 
—  mm 

N  ft 

90 


N 

-El 

i=l 


cos . 


(ei 


V 


•  Tv1- 

“  2N  mln  4 _ 

2N  0- -  •  Li=l 


cos  2(0.  - 


V  +N 


Suppose  is  the  point  whose  angular  coordinate  is  20^  ,  i=l, 

2,  ...  ,  N  .  Let  R  , -of  length  R, .be  the  resultant  of  the  unit 

cos  2(0_^  -  0^)  is  the  component  of 

* 

R  on  A  ,  the  vector  whose  angular  coordinate  is  20-  .  Now  R  is 

~  ~  u 
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fixed  for  a  given . sample; -clearly. . X. -is -minimized  when  A  is  along  the 

direction  opposed  to  .R  . , -so  that  -  X  =r  -rrR  ,  and 

(8)  Js(l  -  R/N)  . 

*Kj 

Similarly  S^x  =  %(1  +.R/N) .  r  when  A  -  lies  along  R  .  If  the  set  of 

points  P  is-uniformly-distributed-over:the:circle,  so  is  the  set  , 

and  R  is  then  the  resultant :  of  ..a  -set.  of :  randomly  distributed  unit  vectors; 

the  distribution  of  its. length  R  . is  well-known  (see;  e,g.;  Greenwood  and 

Durand  (1955) ,  and  Stephens  (1969) . for  critical  values) ,  so  that  critical 

values  of  S  .  and  S  can  be  found  exactly.  Because  of  the  identity 
mm  max 

S  .  +  S  =  1  ,  tests  based  on  these  statistics  will  give  the  same  result, 

mm  max 

as  noted  above.  When  A  is  known,  the  test  statistic  is  > 

S--  =  N  ^  V  cos2  (0 ^  - -0q)  .==  %(H-^X)  ;tcritical/points  for  X/N  are 

given  in  Stephens. (1969) ,  where  the  tests  are  given  in  greater  detail. 

Power  tables  are  also  given. 

If  we  apply  to  the .circle  the. arguments  given  above  for  the  sphere, 
we  find 


(9) 


S 

max 


f 


where  r^  is  the  characteristic  root.of -the  matrix  X  defined  in 
Section  3.1  for  p=2  .  Then . the  density^of  the  limiting  distribution  of 


r^  is  obtained  from  Section .3.1  as 


-4r 


- 8r^  e  . 


2  2 

It  then  follows  that  8r^  has  a  limiting-  X2  ~  distribution.  It  has  also 

2  2 

been  shown  that  for  N  large,  2R./N  .has  the  X2 ~ distribution,  if  we 
write  (8)  as 
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we  see  at  once  the  equivalence, -for. large . .N: of . (8)  and  (9).  However, 
(8)  leads  to  exact -distributions, -and -we: have  no  parallel  exact  results 
for  the  sphere.  It . is -clear -that -since: the  distribution  of  (£,  m,  n) 
does  not -depend -on  the -coordinate -system, :the: exact  distribution  of  the 
entries  of  X  in  (3) -does -not-depend  -on -the: coordinate  system  and  one 
must-be  able-to  write. the-density _in:terms :of  the ; latent  roots,  which 
will  not- change  with  an. orthogonal- transformation  of  coordinates;  if  this 
could  be  found  exactly.  Theorem. 13. 3.1_ of  Anderson  (1958): could  again  be 
used,  ..but  -  the.  problem  appears  -  to -be -a:  difficult  one. 
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Table  1 

LOWER  TAIL  ..PERCENTAGE  POINTS  FOR  S  .  AND  UPPER  TAIL  POINTS 

mm 


FOR  S  ON  SPHERE 
max 


A 

.a  (%)  1.0 

S  .  :  lower  fail 
mm 

2.5  5.0 

10.0 

10,0 

S  :  upper  tail 

max 

5,0  2.5 

1.0 

5 

0.007 

0.011 

0.019 

0.031 

0.714 

0,751 

0.784 

0.821 

6 

.016 

.024 

.034 

.050 

.678 

.712 

.743 

.779 

7 

.026 

.038 

.050 

.067 

,651 

,684 

.712 

.746 

8 

.037 

.051 

.  064 

.081 

.630 

» 661 

.687 

.718 

9 

.048 

.062 

,075 

.093 

.610 

.641 

.667 

.694 

10 

.058 

.073 

.087 

.105 

,596 

.625 

.650 

.677 

12 

0.076 

0.091 

0,106 

0,123 

0,574 

0.598 

0,621 

0.648 

14 

.091 

.107 

,120 

.137 

.554 

.578 

.599 

,623 

16 

.103 

.120 

.133 

.149 

.538 

.559 

.581 

.604 

18 

.114 

.131 

.144 

.158 

,526 

.544 

.566 

.587 

20 

.124 

.140 

.152 

.167 

.515 

.535 

.553 

.575 

25 

.144 

.158 

,170 

.  184 

.496 

.512 

.530 

.550 

30 

0.160 

0.172 

0.184 

0.196 

0,479 

0.495 

0.510 

0.528 

40 

.183 

.192 

o  203 

.214 

.459 

,473 

.487 

.501 

50 

.  .198 

,207 

.216 

.227 

.447 

.460 

.471 

.484 

60 

.208 

.217 

.226 

,235 

,438 

,449 

.458 

.470 

70 

.216 

,226 

,234 

,243 

.429 

.439 

.448 

,461 

80 

.223 

.231 

.239 

.248 

,423 

.432 

.441 

.452 

100 

.233 

.242 

.248 

.257 

,413 

.422 

,430 

.440 

b: 

-1.038 

-0.948 

-0.873 

-0,788 

0.788 

0.873 

0.948 

1.038 

For  N  >  100  ,  use  the.  approximation  1/3  +•  b/xN  ,  where  b  is  given  in  the 
last  row. 


Table  2 

PERCENTAGE  POINTS  FOR.  (MONTE  CARLO) 


Lower  Tail 

Upper 

Tail 

N 

a(%)  1.0 

2,5 

5.0 

10.0 

50.0 

10,0 

5.0 

2.5 

1.0 

5 

0.130 

0.159 

0.183 

0.211 

0.304 

0,390 

0.408 

0.428 

0.444 

10 

.203 

.221 

.238 

.255 

,320 

.379 

,395 

.409 

.422 

20 

.246 

.258 

.269 

.282 

,325 

.368 

.  380 

.389 

.399 

50 

.279 

.  286 

,294 

.302 

.329 

.356 

,364 

.370 

.377 
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Table  3 

POWER  OF  TEST  -FOR. .RANDOMNESS 


The  table  shows  the. proportion -of  Monte  Carlo  samples  significant, 
when  tested. for  randomness  with. test  of  size  a  ,  when  the  sample  of  size 
N  ,  were  in  fact  drawn  from.an. equatorial -distribution  with  given  K  . 


a  = 

0.05 

a  = 

0.01 

K-:  1 

2 

3 

4 

1 

2 

3 

4 

N 

6 

0.060 

0.09 

0.15. 

0.21 . . 

0.01 

0.03 

0.04 

0.05 

10 

j .  072 

.16 

.30.  .. 

.49 

.01 

'  .05 

.09 

.27 

15 

.091 

.27 

.53 

.75 

.03 

.11 

.25 

.49 

20 

>.120 

.40... 

.68 

.92 

.04 

.16 

.40 

.71 

30 

.200 

.58 

.89 

.99 

.05 

.30 

.76 

.95 

50 

.35 

.82 

.99 

1.00 

.11 

.73 

.97 

.99 
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